The QCD Mixed Phase: Shaken, not Stirred 
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Abstract 
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(N 

;_j ■ Baryon fluctuations exceeding Poisson expectations can signal a nearly first order 

phase transition at RHIC. We show how these fluctuations can be measured, and 
apply a dissipative-hydrodynamic formulation used in condensed matter physics to 
00 ■ simulate their evolution. 

(N 

\^ • If the phase transition from quark matter to hadron gas is first order, matter at 

. the appropriate temperatures and densities can exist as a mixed phase consisting 

I of plasma droplets in equilibrium with a surrounding hadronic fluid. If formed in 

' ion collisions, this mixed phase can produce large event-by-event fluctuations as 

the system hadronizes. Extraordinary baryon fluctuations [1] can accompany a first 
order transition at high baryon density [2] and, possibly, a near transition at zero 
I baryon density [3], [4]. In ref. [1], we argued that baryon number conservation and 

rapid longitudinal expansion limits the extent to which post-hadronization interac- 
O , tions can erode fluctuations in a rapidity interval. Here, we further explore the rise 

^ I and fall of these fluctuations using real-time lattice simulations [5] . 

At high baryon density, QCD with two massless flavors can exhibit a first order 
^ ' transition whose coexistance curve culminates in a tricritical point at temperature 

■ Tc and baryon chemical potential fic [2]. For T > Tc and fi < ^ic, ^ second order phase 

transition breaks/restores chiral symmetry. If the quark masses are sufficiently large, 
the second order transition is replaced by a smooth transformation (chiral symmetry 
is explicitly broken). The first order line remains, however, with the tricritical point 
replaced by a critical point in the same universality class as a liquid-gas transition. 

At RHIC, baryon density may also serve as an approximate order parameter for 
the nearly first order transition at small net baryon density. Lattice simulations 
[3] and general arguments [4] show that the baryon susceptibility x at ^ = can 
increase suddenly as temperature is increased near Th ~ 150 MeV, where the chiral 
order parameter and the energy density change sharply. Jumps in the susceptibility 
commonly accompany first order transitions. For a liquid-gas transition, x = dp/dfj, 
is proportional to the compressibility: steam is much more compressible than water. 

Large fluctuations in baryon number occur during phase separation in a first order 
transition. Figure 1 (left) shows the phase diagram in the T — p plane, where p 
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Fig. 1. Phase diagram (left) and free energy (right) vs. baryon density for (3). 



is the baryon density. A uniform system quenched into the outer parabohc region 
below Tc will separate into droplets at high baryon density pq surrounded by matter 
at density p^. The net baryon number Nb in a sub- volume of the system varies 
depending on the number of droplets in the sub- volume. The variance of the baryon 
number V = (A^|) - {Nb)'^ can exceed the equilibrium expectation by an amount 

^V^f{l-f){^NB)\ (1) 

where / is the fraction of the high density phase in the sub-volume V and IS.Nb = 
{Pq ~ Ph)y ■ We will argue that nonequilibrium evolution in ion collisions can allow 
these fluctuations to survive post-hadronization evolution. 

We stress that a super-poissonian variance such as (1) is straightforward to test 
experimentally by measuring 

% = {Vp - (iVp))/(iVp)^ (2) 

where Np is the number of protons in a rapidity interval and Vp is its variance. 
This quantity vanishes in equilibrium and is related to the more familiar scaled 
variance oOp = {Np){l + fip). Significantly, we find that Q.p equals the total VIb 
- which includes unseen neutrons - for a range of thermal and Glauber models 
that respect isospin symmetry. Specifically, Vip = Q,b because the probability of Np 
satisfies p{Np) = "^j^f^ p{Nb)p{Nb\Np) with a binomial distribution p{NB\Np) for 
Np at fixed baryon number Nb- Isospin fluctuations can alter p{NB\Np) near the 
tricritical point or in the presence of a disoriented chiral condensate, but that will 
be evident from pion measurements. 

To describe the evolution of the inhomogeneous mixed phase, we follow the standard 
condensed matter practice and write a Ginzburg-Landau free energy: 

f = K{Vpf/2 + fo, h = -m\p-p,f/2 + X{p-p,f/A (3) 

where /o(/o) describes the excursions of the baryon density p from its equilibrium 
value in the uniform matter. The k term describes the droplet surface tension, 
a oc k}/"^ . For m? (x Tc — T we find the correct liquid-gas critical exponents. The 
values ph and pq in fig. 1 (right) correspond to the equilibrium densities at T < Tc. 
To describe the dynamics of the system, we must account for the fact that baryon 
number is conserved. Furthermore, it is crucial to include dissipation to describe 
this strongly fiuctuating system. The simplest equations that meet these criteria 
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Fig. 2. Order parameter in the transverse plane without (a) and with expansion 
(b). 

are: 

dp/dt = MvV, = /o - '^vV; (4) 

model B in [6]. We identify D = 2m^M as the baryon diffusion coefficient by 
linearizing (4) about ph- 

To describe nuclear collisions, we extend (4) to include drift due to Bjorken longi- 
tudinal flow: 

dp/dr + p/t = MVV, (5) 

where r is the proper time and p is given by (3, 4). The new term forces the average 
density to decrease as (p) oc r~^, driving the system through the phase coexistence 
region. Fluctuations grow when densities are near pc (c.f. fig. 1). 

Phase separation is most dramatic if the rapid expansion drives the system into 
the unstable region; i.e., the inner parabolic region in fig. 1 (left), corresponding 
to fl){p) < 0. Droplets form from runaway density fluctuations in a process known 
as spinodal decomposition. Linearizing near pc, we estimate the time scale for this 
process to be = where ^ = k^I"^ jm is the correlation length. For times t ^ 

the system undergoes a nonlinear evolution in which droplets merge, reducing 
their surface energy. 

Figure 2 compares 2+1 dimensional numerical simulations of (4) and (5) for = 
{p — Pc)/ Pc in the transverse plane. We take tr and ^ to each be 1 fm as motivated in 
[1]. (Z) ~ 8 fm is consistent with calculations in [7]). The expanding system reaches 
Pc at To = 5 fm. Expansion shown in fig. 2b prevents droplets from merging as in 
fig. 2a. Because this is a dissipative system, we must apply thermal noise at each 
lattice site at tq to seed phase separation. The memory of the initial conditions is 
essentially lost for t — tq > tr. 

Figure 3 shows the computed variance for two different initial times and for two ra- 
pidity intervals. The variance is computed from a sample of 5000 simulated events, 
each unique due to the thermal noise. We see that the super-poissonian fluctua- 
tions grow appreciably by r ~ 2tq. This variance drops as the rapidity interval is 
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Fig. 3. Enhanced variance vs. time for two rapidity windows, eq (1). 

increased. We find that variance is governed by the ratio tq/t/j, which compares the 
expansion and droplet-growth time scales. 

We emphasize that these calculations include diffusion, which dampens the fluc- 
tuations once the system becomes stable. For (3) with {p) oc r~^, the system is 
unstable only for r < 2.3 tq. We extend the calculations to much longer times to 
demonstrate that the fluctuations in rapidity survive well past the freezout time, 
of order 10-30 fm, in accord with [1]. We comment that convection, viscosity and 
collision-geometry effects can reduce Op compared to fig. 3. Moreover, our phase 
transition effect may be compensated to some extent by the effect discussed by 
Koch and Asakawa in these proceedings, which owes to the difference between fluc- 
tuations in a plasma compared to a hadron gas. Nevertheless, the strength of the 
signal in our exploratory calculations invites further work. 
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